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Scattering of cosmic microwave background (CMB) radiation in galaxy clusters induces a polariza- 
tion signal proportional to the CMB quadrupole anisotropy at the cluster's location and look-back 
time. A survey of such remote quadrupole measurements provides information about large-scale 
cosmological perturbations. This paper presents a formalism for calculating the correlation func- 
tion of remote quadrupole measurements in spherical harmonic space. The number of independent 
modes probed by both single-redshift and volume-limited surveys is presented, along with the length 
scales probed by these modes. In a remote quadrupole survey sparsely covering a large area of sky, 
the largest-scale modes probe the same length scales as the quadrupole but with much narrower 
Fourier-space window functions. The largest-scale modes are significantly correlated with the local 
CMB, but even when this correlation is projected out the largest remaining modes probe gigaparsec 
scales (comparable to the CMB ai I — 2-10) with narrow window functions. These modes may 
provide insight into the possible anomalies in the large-scale CMB anisotropy. At fixed redshift, the 
data from such a survey form an i5-type spin-2 field on the sphere to a good approximation; the 
near-absence of B modes will provide a valuable check on systematic errors. A survey of only a few 
low-redshift clusters allows an independent reconstruction of the five coefficients of the local CMB 
quadrupole, providing a test for contamination in the WMAP quadrupole. The formalism presented 
here is also useful for analyzing smaller-scale surveys to probe the late integrated Sachs- Wolfe effect 
and hence the properties of dark energy. 

PACS numbers; 98.80.Es,98.70.Vc,95.85.Bh,98.65.Cw 



I. INTRODUCTION 

Cosmologists are making rapid progress in our un- 
derstanding of the structure of the large-scale Universe. 
On the largest scales, the chief source of information is 
the cosmic microwave background (CMB) anisotropy and 
polarization, particularly the all-sky data from WMAP 
[1 H H H lE H' There have been tantalizing hints of 
unexpected behavior in the largest-scale modes of the 
CMB anisotropy The COBE DMR detected a lack of 
anisotropy power on the largest angular scales 0, S 13 ' 
and WMAP has confirmed this result 0, Q, . There is 
evidence suggesting that these largest-scale modes are in- 
consistent with statistically isotropic theories because of 
correlations between modes and /or asymmetry between 
hemispheres [lS|ll|ll|ll|lllil|ll|ll|l3; 

however, 

there is disagreement over how to interpret these results 
0,0,0]. In particular, these results are subject to the 
classic problem of a posteriori interpretation of statistical 
significances: if an unexpected anomaly is found, and its 
statistical significance is computed thereafter, one can- 
not necessarily take the significance at face value. (After 
all, in any large data set, something unlikely is bound to 
occur.) 

The best way to resolve this situation is of course to ob- 
tain a new, independent data set probing the same phys- 
ical scales. Unfortunately, large-angular-scale CMB ob- 
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servations are already at the "cosmic variance" limit, and 
other independent probes of these ultra-large scales are 
few. Observations that may provide independent infor- 
mation on these scales are therefore of considerable inter- 
est. Large-angular-scale CMB polarization data provide 
some relevant information [23, [sj] , although the number 
of independent modes probed is small and the results 
may depend on the details of reionization. 

The scattering of CMB photons in clusters of galaxies 
may shed light on this puzzle. This scattering induces 
a polarization signal 22|, which is determined by the 
quadrupole anisotropy in the photon distribution at the 
cluster location. This "remote quadrupole" signal probes 
large-scale modes of the density perturbation field that 
are different from those probed by the local CMB, so by 
measuring these remote quadrupoles it may be possible to 
get around the cosmic variance limit . It has therefore 
been proposed that a survey of remote quadrupoles may 
shed light on the puzzle of large-scale CMB anomalies 

MM- 

The remote quadrupole signals from different clusters 
are strongly correlated with each other and with the lo- 
cal CMB anisotropy 26] . It is therefore not obvious how 
to design a survey to obtain the maximum amount of 
new information. In addition, we wish to know what 
physical scales of perturbation are probed by a given 
survey; this will depend on both the redshifts and the 
angular distribution of clusters observed. In this article 
I will develop a formalism for determining the indepen- 
dent fluctuation modes that are probed by a survey of 
remote quadrupoles. For a survey that sparsely covers a 
large area of sky, the largest-scale modes probe compa- 
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rable length scales to the first few CMB multipoles, with 
Fourier-space window functions that are narrower than 
that of the local CMB. Determination of these modes 
may be expected to provide insight into the interpreta- 
tion of the possible anomalies in the large-scale CMB 
observations. The value of a sparse large- area survey for 
this purpose has been noted elsewhere [^l. This paper 
provides the first detailed assessment of the amount of 
information available in such a survey. 

The correlation function of remote quadrupole mea- 
surements is quite complicated, depending on both the 
clusters' redshifts and their angular separation [2^. At 
fixed redshift, the remote quadrupole is a spin- 2 field on 
the sky, so it is natural to express it as an expansion in 
spin-2 spherical harmonics. Because it is predominantly 
derived from scalar perturbations, at any given redshift it 
contains (to a good approximation) only E modes, with 
no B contribution. This should provide a valuable check 
on systematic errors in any future survey. 

At low redshift, the measurements naturally become 
strongly correlated with the local CMB temperature 
quadrupole. As a result, the five coefficients a2m of the 
local quadrupole can be easily measured from a survey 
of only a few low-redshift clusters . 

The spherical harmonic basis diagonalizes the angular 
correlations, giving a sequence of correlation functions 
that depend only on redshift. It is much simpler to de- 
termine and count the independent normal modes in the 
spherical harmonic basis rather than in real space. For a 
survey that covers only part of the sky, of course, the in- 
dividual spherical harmonic coefficients will not be mea- 
sured. However, just as in the case of the local CMB 
we can still use the spherical harmonic basis to count 
the number of modes that can be measured, scaling the 
results by the fraction of sky covered. 

On smaller scales, a remote quadrupole survey provides 
insi ght into the growth of structure in the recent past 
[23. l25l |2^ . The remote quadrupole signal, like the 
local CMB, contains contributions both from the surface 
of last scattering and from the integrated Sachs- Wolfe 
(ISW) effect resulting from time variations in the gravi- 
tational potential along the line of sight (See, e.g., 
for an overview of the physics of CMB anisotropy.) 
Since the ISW contribution to the remote quadrupole 
measurements differs from that of the local quadrupole, it 
is possible to extract information about the recent growth 
of perturbations. The formalism developed in this paper 
provides a method of quantifying the amount of extra 
information that can be obtained from such a survey. 

If a cluster has a peculiar velocity, then there is a kine- 
matically induced polarization signal as well as the sig- 



nal considered here [sil . Is^ . This kinematic contribu- 
tion can be removed through multifrequency observations 
|27j . and will be ignored in this paper. In addition, we 
will not consider the polarization induced by scattering 
off of diffuse structure 's^; rather, we will envision a 
survey directed at specific clusters of known redshift. 

This paper is structured as follows. SectionlTTldevelops 
the formalism for calculating the correlation function of 
remote quadrupole measurements. Section IITTI shows the 
information that can be obtained in hypothetical remote- 
quadrupole surveys on a shell at a single redshift as well 
as in volume-limited surveys, and section Hvl contains a 
discussion of the significance of these results. Some more 
than usually boring mathematical steps are contained in 
an appendix. 



II. FORMALISM 
A. Remote quadrupole in a single cluster 

We will assume a flat spatial geometry and label any 
cluster's comoving position with an ordinary 3-vector r, 
with spherical coordinates (r, Oi.,(f>f) defined in some fixed 
earth-centered coordinate system. 

For any particular cluster, we will find it convenient 
to introduce a second coordinate system, denoted by a 
prime, which will have its z' axis is aligned with f, the 
direction from earth to the cluster. To be specific, let the 
primed coordinate system be obtained from the unprimed 
by rotating through an angle Oj. about the x axis and then 
by an angle 0f about the (original) z axis, with the third 
Euler angle set to zero. 

Suppose that an observer in that cluster at the clus- 
ter look-back time measures the CMB anisotropy, conve- 
niently recording the results using the primed coordinate 
system: 

AT 

— {h') = J2aim{r)Yi„,ih'). (1) 

l,m 

The quadrupole spherical harmonic coefficients are 

a2™(r)- y"d3fcA2(fc;r)<5*(k)e^'^ "-r;„(k'). (2) 

Here (5$ is the Fourier-space perturbation in the gravita- 
tional potential. On the large scales of interest to us, the 
quadrupole transfer function contains Sachs- Wolfe and 
ISW terms: 



A2(fc;r) = -y (^j2[M^/-r?rcc)]+6^'' d^'j2[A:(r7-?7')]^(f^))- (3) 
In this expression j2 is a spherical Bessel function, a is the scale factor normalized to unity today, 77 = 770 — r is 
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conformal time (drj = dt/a{t)), rjj-cc is the conformal time of recombination, rjo is conformal time today, and D is 
the matter perturbation growth factor normahzed to unity at high redshift (e.g., |3^). We assume that the dark 
energy is spatially uniform, e.g., a cosmological constant. In addition, we assume instantaneous recombination and 
ignore reionization. Most of the quadrupole signal seen by observers in the cluster is due to photons that come from 
last scattering (just as most of the local quadrupole signal is), so neglecting reionization is a good approximation in 
this context. On the other hand, we cannot ignore reionization when considering the correlation between the remote 
quadrupole signal and the local CMB polarization quadrupole, as discussed below. We work in units where c = 1. 
The observed cluster polarization signal is proportional to the m = ±2 spherical harmonic coefficients: 

p±{v) = {Q±iU) (r) = iVa2±2 (r) . (4) 

where 



and r is the cluster optical depth. We want to study the behavior of p± as a function of cluster position r. Since p± 
are complex conjugates of each other, we need only compute one of them. Let's focus on p_(r), which we will call 
simply p{r) from now on. 
The observed signal is 

p{r) = Na2-2{v) = N ^ d^fc A2(fc; r)5$(k)e^'' '-y22(k')- (6) 
Using equation (|A3p . we can write 1^22 (k') in the unprimed coordinate system: 

p{v)=nJ^ J2 (-l)™2r2_™(f) / d'kA2{k;r)S^{k)e'''-^Y2M, (7) 



where 25^2-m is a spin-2 spherical harmonic. 

For a fixed distance r, p is a spin-2 function of direction f, so it is natural to expand in spin-2 spherical harmonics: 



p{r) = J2PLM{r) 2YLM{r), (8) 




with coefficients given by 

PLAiir) = J dHpiv) 2YIm (9) 

d3A:A2(fc;r)<5*(k)y2m(k) j dH 2Y2-m{v) 2F£M(f )e*'^ "". (10) 

By expanding the exponential in spherical harmonics as shown in the Appendix, we can express the coefficients in 
the following form: 

PLM{r) = i"^ J d''kA2{k;r)6^{k)FL{kr)Yl^j{k), (11) 

where 

FUx) ^V20^ (-1)^'~"^/'(2A + 1) (l ^2 o) (o o) (^2) 

\=L-2,L,L+2 ^ ^ ^ ^ 

It is straightforward to check from equation Ullfl that all coefficients pim{r) are real. In the terminology of CMB 
polarization [sslls^ . this means that the remote quadrupole data form an i?-type spin-2 field at any given distance, 
with no B modes. The absence of B modes arises because we have considered only scalar perturbations as the 
source of the CMB quadrupole. If tensor perturbations were included, then in principle a B component would arise. 
Considering the difficulty of detecting a remote quadrupole signal at all, the prospect of searching for a subdominant 
-B-type signal sounds extremely daunting. It is probably more realistic, therefore, to search for the scalar (£^-type) 
signal in such a survey, using the predicted absence of B modes as a check on systematic errors and noise (see Section 
Ml- I 

B. Correlations between clusters information that can be obtained from such a survey, we 



Suppose that many clusters have been observed at 
many different positions r^. To determine the amount of 
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FIG. 1: Correlation functions for I = 2,3,4,20. 



need to know the correlations {p{ri)p{rj)) . The full cor- 
relation function depends on the directions as well as the 
distances of the clusters |2^. The correlation functions 
are simpler to work with in spherical harmonic space, as 
the orthogonality of the spherical harmonics implies that 
there are no correlations between plm and pl'M'- 

{PLM{r)phM'ir')) ^ £,L{r,r')6LL'SMM'- (13) 
The correlation function is 

Ur.r')^{pLMir)plM{r'))^ 

/ dkk^A2{k;r)A2{k;r')FL{kr)FL{kr')P^{k), (14) 
Jo 

where the power spectrum is given by 

(<5$(k)<5;(k')) = F*(fc)<5(3)(k - k'). (15) 

On the scales of interest, P$(fc) c>c where n ~ 1 is 

the spectral index. 

Fig. n shows the correlation functions £,i for several 
values of In calculating the ISW integral and in con- 
verting from distance r to redshift z, a flat Friedmann- 
Robertson- Walker cosmology with 51™ = 0.3 and JIa = 
0.7 was assumed. The results are normalized to WMAP 
on large angular scales. At low I, the correlations are 



extremely broad, meaning that even a survey covering 
a wide range of distances will contain few independent 
modes per {l,m), as we will see in section Fill Bl 

Remote quadrupole measurements are in general cor- 
related with the local CMB anisotropy and polarization. 
To assess how much independent information is con- 
tained in the remote measurements, we need to know 
how strong these correlations are. Let us begin by con- 
sidering correlations with the locally-measured temper- 
ature anisotropy. If the anisotropy spherical harmonic 
coefficients are aim, then the cross-correlation is 

air) ^ {PLM{r)aU) (16) 




dk fc2A2(fc; r)AL{k; 0)FL{kr)P^{k), (17) 



where is the appropriate transfer function. 

If we wish to study only the information in a remote 
quadrupole data set that is independent of the local 
anisotropy, then we should project the mode coefficients 
onto the space orthogonal to that probed by the local 
CMB: 

PLAiir) = PLAiir) - ^l^flLM, (18) 

where Ci = {aima*im) is the usual CMB angular power 
spectrum. 
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FIG. 2: Correlation functions for 1 — 2. In the left panel, only the correlation with the local temperature quadrupole has 
been projected out. In the right panel, the correlation with the polarization quadrupole has also been removed. 



We can define a correlation function (^f and perform 
a similar projection to remove the portion of the signal 
that is correlated with the CMB polarization anisotropy 
coefficients af^y^- (There is no significant correlation with 
the i?-type CMB polarization.) In performing this pro- 
jection, it is important to include the effects of reioniza- 
tion, as most of the low-Z polarization signal comes from 
post-reionization scattering. 

Fig.[21shows the correlation functions corresponding 
to the projected coefficients for 1 = 2. As we will see, 
the difference between and £^j^ decreases at higher Vs. 
In performing the projection, the Universe was assumed 
to have completely reionized at z = 11, but the results do 
not depend strongly on the details of reionization, unless 
there was considerable patchiness on large scales. 

The r.m.s. power due to all modes at a given / is 



21 + 1 
An 



(19) 



The real-space correlation functions are easily com- 
puted from the spherical harmonic space functions. The 
correlation between remote quadrupole signals of two 
clusters at locations ri , r2 is 

°° 2r + 1 

(p(ri)p(r2)) = a(ri,r2)PL(fi •£•2) (21) 



L=2 



using equations (|HJ| and Ijl^fl and the spherical harmonic 
addition theorem. Similarly, the correlation between a re- 
mote quadrupole measurement p(ri) and the local CMB 
(AT/T)(r2) is 

°° or I 1 

(p(ri)(Ar/T)(r2)) = ^ —^CL{n)PL{h-h). (22) 



L=2 



Once the initial investment of calculating the /-space cor- 
relation functions has been made, these formulae allow 
rapid calculation of real-space correlations. 



and similarly for p^. As Fig. |3| indicates, the quadrupole 
I = 2 dominates the unprojected power at low redshift, 
but power shifts to smaller angular scales (higher I) as 
the redshift increases. At / = 2, the projected power 
is much less than the unprojected power at all redshifts: 
the modes p2m are strongly correlated with the local tem- 
perature quadrupole at low z and with the polarization 
quadrupole at high z. Modes with I > 2 are compara- 
tively weakly correlated with the local signals over some 
redshift ranges. 

All correlation functions except ^2 go to zero as 
z, z' 0. This is expected: at low redshift, the remote 
quadrupole p{r) contains precisely the same information 
as the local quadrupole coefficients a2m , so it must trans- 
form as a quadrupole itself. Indeed, it is straightforward 
to check from equation (|1HI that 

Pirn{r) N \ —aim5i2 as r ^ 0. (20) 



III. SCALES PROBED BY REMOTE 
QUADRUPOLE SURVEYS 

A. Survey at a fixed redshift 

We next examine the length scales probed by the vari- 
ous multipoles, assuming an all-sky survey has been used 
to estimate the coefficients pim at some fixed redshift. We 
can write the signal as 

pf{r) = J dkk^P^{k)W^{k;r), (23) 

with Wi{k;r) = ^A2{k;r)Fi{kr). Since cx k~^, 
the quantity W^/k is proportional to the power per 
wavenumber interval dk. Similarly, we can define a win- 
dow function for the quantity p/" that results from pro- 
jecting out the part of the signal that is correlated with 
the local CMB. 



6 




FIG. 3: The r.m.s. signal per multipole. The left panel shows 
the signal in modes / = 2 (solid, black), 3 (dashed, red), and 
4 (dot-dashed, green) as a function of redshift. The bumps 
at < 0.5 are due to the ISW effect. In each case, the three 
curves from top to bottom are the total signal, the result of 
projecting out the correlation with temperature quadrupole, 
and the result of projecting out both temperature and polar- 
ization. The right panel shows the signal as a function of / for 
redshifts z = 0.1 (solid, black) 1 (dashed, red), and 2 (dot- 
dashed, green). The upper curve is total signal, and the lower 
curve is the result of projecting out the correlations with local 
temperature and polarization. 



The first few window functions are siiown in Fig. ^ 
for z = 2. Window functions corresponding to both pi 
and are shown. The range of scales probed by the 
various window functions are indicated with horizontal 
error bars, and for comparison the ranges corresponding 
to the local CMB power spectrum C2 and Cio are also 
indicated. Because of the ISW effect, the local CMB 
window functions are quite broad. 

The first few unprojected modes probe scales as large 
as the CMB quadrupole but with narrower window func- 
tions. As noted earlier, these are significantly correlated 
with the local CMB polarization multipoles. Nonethe- 
less, considering the likelihood that large-angle CMB 
polarization multipoles may be contaminated by fore- 
grounds or systematic errors, the unprojected modes will 
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FIG. 4; Window functions. The quantity /k, which is 
proportional to the mean-square power in multipole a given 
multipole, is shown for a survey at z — 2 for 1 = 2 (solid, 
black), 3 (dashed, blue), and 4 (dot-dashed, red). In each 
case, the left (thicker) curve is the window function for the 
total signal pim, and the right (thinner) curve is the window 
function for p^. All functions are normalized to integrate to 
one. The bars above each window function show the 25th, 
50th, and 75th percentiles. The green bars at the top of 
the plot indicate the percentiles for the I = 2 and I = 10 
multipoles of the local CMB anisotropy. 

still provide valuable new information on the largest-scale 
perturbations in the Universe, or at least test our under- 
standing of large-angle polarization data. (The remote 
quadrupole survey will of course be susceptible to sys- 
tematic errors and foregrounds as well, but the suscepti- 
bility will be different from that of the local polarization 
multipoles.) 

The projected modes probe smaller scales, but they 
are still in the gigaparsec range, comparable to the first 
10 or so CMB multipoles, and in some cases have narrow 
window functions. In practice, the projected modes at 
I = 2 {P2) are unlikely to be reliably measured, because 
the correlations are so strong, but projected modes with 
I > 3 will allow us to probe these large scales. 

Fig.|Slshows the r.m.s. power pi and per multipole, 
plotted against the effective scale for each multipole. In 
interpreting this plot, bear in mind that each point rep- 
resents the r.m.s. signal from all 21 + 1 modes at a given 
I. 

In order to measure the quantities p^, in principle 
we need an all-sky cluster survey, knowledge of the lo- 
cal CMB anisotropy and polarization spherical harmonic 
coefficients, and knowledge of the correlation functions 
Q{r) and C^{r) in order to project out the local con- 
tribution. In practice, of course, difficulties are likely 
with all of these. Section Hvl contains some discussion of 
how to mitigate these problems. For the moment, ob- 
serve that information on large physical scales is found 
at large angular scales. We must survey a large fraction 
of the sky if we want to address the puzzles in the large- 
scale CMB with this technique. However, note that at 
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FIG. 5: The r.m.s. signal vs. effective length scale. The upper plots are for a survey at 2 = 2, and the lower plots are for z = 3. 
On the left the total signal pi is shown, and on the right is the projected signal p/". The horizontal error bars indicate 25th, 
50th, 75th percentile contributions to the signal. Multipole I increases from left to right within each plot. The dashed green 
bars at the bottom indicate the window functions for the local CMB anisotropy at Z = 2 and I — 10; the vertical position of 
these bars is arbitrary. 



redshifts z — 2-3 the signal drops fairly rapidly as a func- 
tion of I. This is good news: it means that a relatively 
sparse survey can measure the low-/ modes without ex- 
cessive contamination from small angular scales. 



B. Volume-limited survey 

In the previous subsection, we considered surveys at a 
fixed redshift. We now imagine a volume- limited survey 
out to some maximum redshift Zmax- Let us continue to 
assume an all-sky survey, so that it is natural to think of 
the survey in spherical harmonic space. In this case, our 
survey provides estimates of each of the functions pim (r) 
at multiple values of r. 

For each I, we can enumerate a list of signal strength 
eigenmodes tpni that are solutions to 




^iir,r')^bniir')r'^dr' = Xni^r). (24) 



The mode functions ipniir) 2^m(f) form an orthonormal 
basis, which we can use to express the signal p{r). The 
mean-square signal in each mode is the eigenvalue A„i , so 
these modes provide a useful guide to tell us where the 
signal is strong. 

For each mode, we can calculate a window function and 
hence assign a range of wavenumbers probed as we did 



for the surveys at fixed redshift. Results are illustrated 
in Fig. El 



As one would expect, the mode with highest signal at 
each I corresponds to a simple weighted average of pim as 
a function of r with positive weight everywhere. The next 
mode is essentially a difference between low- and high- 
redshift signals, and subsequent modes contain more ra- 
dial oscillations. As Fig.|n|indicates, only the first couple 
of modes are likely to be measurable at any given m). 
Once again, the I = 2 modes are strongly correlated with 
the local CMB polarization. Assuming the large-angle 
CMB polarization has been well measured, they provide 
relatively little new information; modes with I > 3 are 
the richest source of independent data on large-scale per- 
turbations. On the other hand, if the first few CMB 
polarization multipoles are uncertain due to foregrounds 
or systematic errors, then the I = 2 modes of a remote 
quadrupole survey may help to fill in this gap. 



A comparison of Figs. |S1 and El shows that signifi- 
cantly more large-scale information can be obtained from 
a volume-limited survey than from a survey on a shell. 
There is, of course, an obvious price to pay: many more 
clusters must be observed to estimate all these modes. 
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FIG. 6: The r.m.s. signal vs. effective length scale. The upper panels show the results of a survey out to a maximum redshift 
of 2, and the lower panels are for a maximum redshift of 3. The total signal pi is plotted on the left, and the projected signal 
p'^ is on the right. The horizontal axis indicates the length scales of each mode as in Fig. |S] Solid black bars are the highest 
signal-strength eigenmode as a function of I. The dashed red bars are the second mode for each I, and the dot-dashed blue 
bars are the third mode. Multipole I increases from left to right within each category. The long-dashed green bars at the top 
indicate the window functions for the local CMB anisotropy at i = 2 and I = 10; the vertical position of these bars is arbitrary. 



IV. DISCUSSION 

We have seen that an all-sky survey can probe the 
gigaparsec-scale Universe, measuring fluctuation modes 
that are independent of the local CMB. In surveys at red- 
shifts around 2-3, the large- angular scale modes provide 
data on perturbations on the same length scales as the 
first few CMB multipoles, but with quite narrow window 
functions. 

The results shown in the previous section were for 
an idealized survey: in addition to full sky coverage, 
the local CMB anisotropy coefficients aim and the cross- 
correlations Q (r) , as well as the corresponding quantities 
for polarization, were assumed to be known in order to 
compute the projected signal p^^^. We must ask what 
happens if these assumptions are replaced by more re- 
alistic ones. The most complete way to answer these 
questions would be to assume a precise survey geometry 
and compute the resulting Fisher matrix. We will not 
perform such a detailed analysis here; we can, however, 
make some general observations. 

In a survey that covers a fraction of the sky /sky, 

— 1/2 

only band powers with width Al ~ f^^^' can be recov- 
ered, not individual multipoles. Furthermore, the lowest- 
l modes cannot be recovered at all. For the goal of prob- 
ing the largest scales, therefore, large sky coverage is es- 
sential independent of the choice of redshift. A survey 



with /sky = 0.1, for instance (4000 square degrees) would 
be able to recover only a single mode in the / = 2-3 band. 

On the other hand, the power drops fairly rapidly as a 
function of Z, so contamination of the \ow-l modes from 
high-/ power is modest. In other words, in order to probe 
large scales, we should survey as much sky as possible, 
but the survey can be sparse. 

Next, let us consider uncertainties in projecting out 
the local CMB contribution (i.e., going from pim to p^)- 
For all I > 2, this projection is subdominant to the pri- 
mary signal over some range of redshifts, so independent 
information should be obtainable from these modes. 

The I = 2 modes are a different matter, as the corre- 
lations are extremely strong there. As Fig. O indicates, 
the projected coefficients are much smaller than the 
unprojected coefficients P2m at all redshifts. At low red- 
shift, the culprit is the temperature quadrupole, while at 
high redshift p2m becomes very strongly correlated with 
the polarization quadrupole. To put the situation pes- 
simistically, accurate extraction of p^„j(r') may never be 
feasible. To extract information about large-scale per- 
turbations that is independent of the local CMB, we will 
look to modes pim with I > 3 (or, in the case of a partial- 
sky survey, by modes that cover the largest available an- 
gular scales but are orthogonal to the quadrupole). 

A more optimistic interpretation is that measurement 
of P2m at a couple of different redshifts can allow us to 
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determine the local CMB temperature and polarization 
quadrupole coefficients (that is, the 5 coefficients a2m 
and the 5 coefficients a^). Since direct measurements 
of these coefficients may be contaminated by foregrounds 
or systematic errors (especially in the case of polariza- 
tion), such an independent determination of these coef- 
ficients will be important in assessing the significance of 
the large-scale anomalies in the CMB. Furthermore, by 
measuring p2m as a function of redshift, we may be able 
to test our theoretical predictions of the cross-correlation 
functions C2, Cl'' ^^vlS providing a probe of the recent ISW 
effect. 

A common question is whether a remote quadrupole 
survey can "beat cosmic variance." The answer depends 
on precisely what we mean by this phrase. Fig. [7| pro- 
vides one possible answer. The figure shows the cumula- 
tive number of independent modes on scales larger than 
a given value for ideal all-sky surveys out to a speci- 
fied redshift. All of the signal eigenmodes are included 
in this count. The number of modes contained in the 
all-sky CMB temperature anisotropy data (without po- 
larization) is shown for comparison. Although a survey 
that went all the way out to 2: = 00 would "beat" the lo- 
cal CMB, realistic surveys never do. The amount of new 
information can be comparable, however, on some scales. 
In particular, the number of new modes obtainable by a 
remote quadrupole survey in the range k ^ 0.002ft,/Mpc 
is about the same as that contained in the local CMB 
(because the slopes of the cumulative curves in Fig.[7|are 
about the same there). Considering the unsettled state 
of our understanding of gigaparsec-scale perturbations 
and the hints that something surprising may be going on 
there, it is clear that there is valuable information to be 
gained. 

This article has focused primarily on the largest-scale 
information contained in remote cluster surveys. The for- 
malism described here is also useful for surveys designed 
to probe the ISW effect HEllllllil . Such surveys pro- 
vide a powerful probe of the recent growth of structure 
and hence may shed light on the nature of dark energy 
and the growth of structure. Because the ISW effect is 
most important at low redshift (see Fig. ^ , such a sur- 
vey will be quite different from those considered here: the 
best approach appears to be a denser survey of a smaller 
area of the sky at low redshift. In planning a survey to 
probe the ISW effect, it will be important to quantify 
the number of independent modes that can be probed. 
The detailed answer will depend on the precise locations 
of the clusters to be surveyed, but a simple estimate ob- 
tained by counting modes in spherical harmonic space 

— 1/2 

and scaling by f^^^ should provide valuable guidance. 

Clusters are of course not randomly distributed "test 
particles" : they are overdensities. One might worry that 
this would lead to biases in the modes recovered from 
such a survey. A remote quadrupole survey (even a small- 
scale one optimized for characterizing the ISW) primarily 
probes scales of several hundred Mpc or more, which is 
considerably larger than the scale associated with the for- 
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FIG. 7: The cumulative number of modes as a function of 
scale. From bottom to top, the solid curves show predictions 
for surveys with Zmax ~ 1, 2, 3. The red dashed curve shows 
the number of modes probed by the local CMB anisotropy, 
and the green dotted curve shows the results of a hypothetical 
remote quadrupole all-sky survey out to 2; = 00. 



mation of individual clusters. One would therefore not 
expect significant bias due to the locations of individual 
clusters. On the other hand, the modes recovered from 
such a survey would presumably be correlated with trac- 
ers of large-scale structure on hundred-Mpc scales. In 
analyzing the results of such a future survey, one would 
want to characterize those correlations, presumably via 
iV-body simulations. For the gigaparsec-scale surveys 
that are the primary focus of this paper, of course, clus- 
ters can be taken as randomly-distributed test particles. 

In a detailed Fisher-matrix analysis of a potential sur- 
vey, the real-space covariance matrix {p{vi)p{Yj)) will be 
needed, as will the correlation with the local CMB. The 
formalism in this paper provides a useful way to com- 
pute these quantities. The full covariance matrix can be 
computed with and without the ISW effect, and Fisher 
matrix estimates of the errors with which ISW parame- 
ters can be reconstructed can be computed quickly and 
easily for any desired survey geometry. 

A survey of the sort considered here will surely be 
a daunting task. The signals are sub-microkelvin, and 
there is unfortunately no shortage of confusing signals. 
Some signals (diffuse Galactic foregrounds and the kine- 
matic signal due to the cluster's peculiar velocity) can 
be distinguished by their spectral signature, assuming a 
multifrequency survey, but detailed information on the 
spectral and spatial properties of polarized foregrounds 
will be necessary. The three- year WMAP data have ad- 
vanced the state of knowledge in this area considerably 
0, and further information will come from the Planck 
satellite Is^ as well as ground-based experiments. 

Of potentially greater concern is the intrinsic CMB po- 
larization (both due to last scattering and reionization) , 
which will be lensed by the cluster itself. In order for the 
remote quadrupole survey to be detectable, we will prob- 
ably need detailed knowledge of the cluster optical depth 
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as a function of position on the sky, and possibly the 
projected mass density as welL With this information, 
we can construct a template for the remote quadrupole 
signal and use it to fit for the two parameters that de- 
termine the remote quadrupole at that cluster [Q and 
U or equivalently the real and imaginary parts of p(r)]. 
Since the background polarization is not expected to be 
spatially correlated with this template, this will help in 
separating the signal from the confusing background. In 
the near future, Sunyaev-Zel'dovich surveys will provide 
a wealth of detailed cluster data [3^ [S^ , so there is rea- 
son to hope that such an approach may soon be feasible. 

In any remote quadrupole survey, assessment of the er- 
rors will be crucial. For instance, errors in determining 
the optical depths t of the clusters can induce spurious 
signals. Just as in the case of CMB polarization maps, 
a valuable diagnostic can be obtained by considerin g th e 
decomposition of the data into E and B modes [S^, ISg • 
At any fixed redshift z, the remote quadrupole data con- 
sist of a spin-2 field on the sphere that is derived from 
a scalar perturbation (assuming that primordial tensor 
perturbations can be neglected). As noted in Section ITTl 
the true signal — everything calculated in this paper — 
should therefore consist only of E modes, precisely as in 
the case of scalar perturbations in the CMB. Noise and 
systematic errors, on the other hand, are likely to popu- 
late E and B equally 'id\. When analyzing results of an 
actual survey, therefore, the B modes can be monitored 
to determine the errors. In practice, for a partial-sky sur- 
vey with sparse sampling, there will be significant E-B 
mixing (41., ,42^ ,433 ,44], but this technique should still 
provide a valuable check. 
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bols. 



1. Rotation matrices and spin-2 spherical 
harmonics. 



For a given cluster location r, we adopt a primed co- 
ordinate system obtained by rotating the z axis until it 
points in the direction f. Let R be the rotation that 
relates the two coordinate systems. The Euler angles as- 
sociated with this rotation are (0f., 6^, 0), using the same 
conventions as '41]. The spherical harmonic Y22 in the 
primed coordinate system can be expressed in the un- 
primed coordinates as 



(Al) 



where D\^^,{R) is the Wigner matrix for the rotation R. 
The Wigner matrices can be expressed in terms of spin-s 
spherical harmonics |4l|: 



The result is 



47r 
21 + 1 



sYi„,{0,^)e''^. (A2) 



Y22{-k') = J —J2h^r2Y2-Ur)Y2M- (A3) 



2. Integrals over spherical harmonics. 



The derivation in Section ITU contains an integral 



APPENDIX A 



/ = Jd^r 2^2-™ (f) 2>7:M(f)e'"'-"". (A4) 



In this section we derive some identities involving 
spherical harmonics, coordinate transforms, and 3-j sym- 



To evaluate this integral, we expand the exponential in 
spherical harmonics to get 



J 



/ = 4^5]z^j;,(fcr)rv(k) / d^fr;^(f) 2Y2-Ui) 2YlMir)- 

A, /J, 



(A5) 



Next, we want to evaluate the integral over the three spherical harmonics. Using the identity sYj^ ~ (•_2)ni+s _^Yi^jn, 
the integral we need can be written in the form 



J= / d^r 2Yi,rnA^) -2Yi^m.,{r)Yi,m,{r). 



(A6) 
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Equation (B3) of tells us how to express all of the spherical harmonics in terms of D-matrices: 



(2;i + l)(2^, + l)(2/3 + l) f , ,3 

Here the ZJ-matrices can be evaluated for any rotation R{4>, 0, 4') with the first two Euler angles being the spherical 
coordinates of f. Since the integrand doesn't depend on the third Euler angle ip, we can replace J cPr with J cfiR, 
an integral over the entire rotation group. Zare |45| (p. 103) gives this integral in terms of 3j symbols: 



(2/i + l)(2Z2 + l)(2;3 + l) ^ h h h \fh h h 



V 47r \-mi -m2 J \2 -2 Q J ' 

So we can write the integral in equation ljA4p as 

^ = E ^'(-l)™JA(fcr)rv(k) V20^(2A + 1)(2L + 1) (^^ ^ . (A9) 

We can use this result to write equation (|10|l as 

PLM{r) = AtiN I d^k A2(fc; r)(5$(k) ^ i V(2i + 1)(2A + l)jx{kr) (^^ K{k), (AlO) 

where 

^« = E Q M i) ^A,(k)r2™(k). (All) 

Expand K in spherical harmonics: ^(k) — X^/omo -^/o^^o^^o^o(^)- "^^^ coefficients are 

= E fm i/ / d^kY^MYx^^imimoi^) (A12) 



/TV"o\-f2 L ,/ 5(2/o + l)(2A + l) ^2 /p X\ f 2 /„ A\ 

^ Z^lm Af J V Att lo oil™ -mo J ^ ^ 



_ (-1)"° 5(2/o + l)(2A+l) (2 lo X\ 

using equation (3.119) in to integrate the product of three spherical harmonics, and then using the orthogonality 
of the 3-j symbols [equation (2.32) in |45|]. So K has only one term in its spherical harmonic expansion: -R'(k) — 
Kl-mYl-mO^)- Substituting this into equation IjAlOp . we get 

PiM(r) = ^^V20^^^\2A + l) (^^ (o o) / d'k A2ik;r)S^{k)Mkr)YlM. (A15) 



The 3-j symbols vanish whenever the triangle inequality Furthermore, ( q o o ) ^ *-* when 2 + L + A is odd. So the 
is not satisfied, so A must be between L ~ 2 and L + 2. sum above contains only three terms: A = L — 2,L,L + 2. 
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